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Abstract 

A numerical investigation for the stability of the incompressible slip flow of normal quan- 
tum fluids (above the critical phase transition temperature) inside a microslab where surface 
acoustic waves propagate along the walls is presented. Governing equations and associated 
slip velocity and wavy interface boundary conditions for the flow of normal fluids confined 
between elastic wavy interfaces are obtained. The numerical approach is an extension (with 
a complex matrix pre-conditioning) of the spectral method. We found that the critical 
Reynolds number {Recr or the critical velocity) decreases significantly once the slip velocity 
and wavy interface effects are present and the latter is dominated {Rccr mainly depends on 
the wavy interfaces). 

PACS numbers : 68.03.-g, 67.57.Np, 67.57.Hi, 67.55.Hc, 67.40.Vs, 43.35.Lq 

1 Introduction 

Ultrasound measurements play very important roles in a '^He investigation^. The time scale of 
ultrasound corresponds to quasiparticle hfe times at low temperatures and first to zero sound 
crossover was observed. Energy scales of ultrasound match the binding energies of Cooper 
pairs in superfluid ^He were observed coupling through density oscillations. Sound transmis- 
sion method using many types of surface acoustic wave (SAW) sensors have been developed 
for acoustical and electrical property measurements of adjacent liquid or gas (especially liquid 
"^He at low temperature^). A Rayleigh-SAW propagates along the substrate surface by emitting 
compressional waves into the quantum fluid and thus the sampling of the Rayleigh-SAW is de- 
termined by the acoustic impedance of the surrounding quantum fluids. Relevant problems are 
the instability for the flow induced by the SAW'^'^ or the slip effect existing along the interface of 
the quantum fluid. In fact, the instability may triggers the laminar-turbulent transition for flows 
of the superfluid. The important issue is to determine the critical velocity^'^ for the relevant 
physical parameters so that the experimental measurements could be under control. 
A most striking characteristic of liquid helium is that it exists in the liquid state down to the 
absolute zero temperature because (i) the van der Waals forces in helium are weak; (ii) the zero- 
point energy, due to the light mass, is large. In fact, it is called a quantum liquid due to these 
kinds of quantum effects and is closely related to the Bose-Einstein condensation for He II. The 
well-known properties of He II can be largely accounted for on the basis of phenomenological 
two- fluid theory^'^. One of the basic assumptions of the two- fluid model is : He II consists of 
a kind of mixture of two components, a normal component and a superfluid component. The 
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former has viscosity while the superfluid can move without friction as long as certain velocity 
limits are not exceeded. 

One crucial issue for the related researches about He II is the critical velocity^'*^'^ (when it 
flows through a capillary or plane channel) which depends on the micro-channel size. Landau 
attributed the existence of a critical velocity in He II to the breakdown of the superfluid due to 
the creation of excitations (he proposed that phonons and rotons axe two types of excitations 
which make up the normal fluid). If the velocity is less than the critical value, there will be 
no dissipation or friction along the boundary or interface for the flow. In plane channel or slab 
flows of quantum fluids, to determine the critical velocity corresponds to finding out the critical 
Reynolds number^~^'^°. 

The traditional starting point of an investigation of hydrodynamic stability is eigenvalue anal- 
ysis, which proceeds in : (i) linearize about the laminar solution and then (ii) look for unstable 
eigenvalues of the linearized problem. In much of the literature of hydrodynamic stability at- 
tention has been restricted to 2D perturbations, in particular, the well-known Orr-Sommerfeld 
equation is an eigenvalue equation for 2D perturbations. For pressure-driven flow, with the 
no-slip boundary conditions (which are normally valid in macro-channels), eigenvalue analysis^^ 
predicts a critical Reynolds number Rscr = 5772 at which instability should first occur, but in 
the laboratory^^ , transition to turbulence is observed at Rscr as low as ^ 1000. 
As for quantum fluids or liquids, there always exists a slip velocity along the static interface or 
confinement due to the microscopically incomplete momentum exchange therein^^. The distinc- 
tion among various flow regimes can be obtained by introducing the Knudsen number (which 
also characterizes the value of the slip velocitylS), Kn=mfp/L, where mfp is the mean free path 
and L is the characteristic flow dimension. Slip flow conditions exist for 0.001 <Kn< 0.1, where 
the flow can be considered continuous if the slip velocity at the walls are taken into account. 
Meanwhile, the effect of elastic or deformable interfaces, like surface acoustic waves (SAW) 
interacted with volume and surface phonons (propagating along the elastic boundaries^^"-'^^), 
upon the stability of quantum fluid flows, however, were seldom considered before, to the best 
knowledge of the authors. Although the mathematical difficulty is essential therein. Note that, 
entrainment of fluids induced by SAW propagating along deformable boundaries have been 
studied-^^ since early 1950s. The role of elastic macroscopic walls resembles that of microscopic 
phonons^. As presented in [4], for the description of the transport processes in the nonequilib- 
rium gas-solid system including the processes occurring in the case of propagation of sound in a 
solid, we need to solve 
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and the associated boundary conditions (scattering and interacting laws near the interface) 



/, n, and H denote the distribution function for gases, volume phonons, and surface phonons, re- 
spectively. Ig, It,, and Is are the corresponding collision integrals. Please see the details in [4] for 
other notations or symbols. The stability problems for such a system of differential-integral equa- 
tions remain open up to now. To escape from above (many-body problems) difficulties, we plan 
to use the macroscopic approach which is a complicated extension of previous approachcs"'^'''"'^^. 
In this work, the stability analysis of incompressible slip flows for a normal fluid (above the tran- 
sition temperature Tc) inside a microchannel with the flow confined by two elastic (wall) layers 
of thin films where surface (acoustic) waves are propagating along the interfaces is conducted. 
Namely, we will relax the static- or rigid-interface boundary conditions which are frequently used 
in previous theoretical and/or experiment approaches into the dynamic- or deformable-interface 
boundary conditions which are more realistic in essence when we consider the flow stability 
problem in a microdomain where surface acoustic waves^^ propagating along the boundaries of 
quantum fluids^'^. The verified code which was based on the spectral method developed by 
Chu^^ will be extended here to include the boundary conditions coming from SAW existing 
along the interfaces to obtain the stability characteristics of the basic flow. 
This paper is organized as follows. We introduce the mathematical formulation for the Orr- 
Sommerfeld equation and the relevant linear stability analysis in Section 2. Boundary condi- 
tions for the slip flow and the wavy interface of the fluid system are then incorporated. The 
expression for primary slip flows of quantum fluids we shall investigate their stability charac- 
teristic will be derived before we describe the numerical approach : a modified spectral method 
at the end of this Section. Two physical parameters, Kq (relevant to the SAW effect) and Kn 
will be introduced or defined. Finally we shall present our calculations and discuss them. Our 
results illustrate that the critical Reynolds number (Rccr) decreases (to 1441) rapidly once the 
interfaces are subjected to propagating waves (or noises) and there are slip velocities existing 
along the boundaries or interfaces {Kq = 1 and Kn=0.001). However, the slip velocity (adjusted 
by Kn) effect is minor and adverse compared to the SAW effect (tuned by Kq) considering the 
decrease of Rccr- 

2 Governing Equations 

Macroscopically, the motion of the normal fluid (above the critical phase transition temperature 
of the quantum fluids) as a whole could be treated by using hydrodynamical models starting 
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from the microscopic atomic wave function^' Here, after the simphfying treatment of the com- 
phcated mathematical derivations, the dimensionless equations of motion for an incompressible 
normal fluid flow^'^'^'^^'^®, in the absence of body forces and moments, reduce to 

^ + (U-V)U = -VP+-^V2U (1) 

where U and P stand for the velocity and pressure of fluids. Re = pumaxh/n is the Reynolds 
number with p and fj, being the density and viscosity of fluids. For the case of normal fluid flows 
driven by a constant pressure-gradient, i.e., plane Poiseuille flow, the length scale is the half 
width of the normal-fluid layer h, and the velocity is the center-line velocity Umax- Following the 
usual assumptions of linearized stability theory^'^"^^, Ui{xi,t) = Ui{xi) + u'i{xi,t), and similarly, 
P{xi,t) = p{xi) +p'{xi,t), the linearized equation, which governs the disturbances are: 

_i + (u . + (u' • V)n, = -Vp' + —VM (2) 

Disregarding the lateral disturbances, w'=0, a stream function for the disturbance, i/j, can be 
defined such that u' = dip/dy, v' = —dil^/dx. Using normal mode decomposition analysis, ip may 
be assumed to have the form ip{x, y, t) = (f){y) ex.p[ia{x — Ct)], a is the wave number (real), C is 
Cr + iCi- This is a kind of Tollmien-Schlichting transversal waves, is the ratio between the 
velocity of propagation of the wave of perturbation and the characteristic velocity, Cj is called 
the amplification factor, and a equals to 27rA~^, where A is the wave length of the Tollmien- 
Schlichting perturbation^^. Substituting the stream function and eliminating the pressure, we 
have the linearized disturbance equation 

_ a2)(L>2 _ ^2-)^ ^ iaRe[{u - C){D^ - a^)^ - {D'^u)(t>] (3) 

where D = d/dy. This is also valid for the slip flow regime^^, 0.001 < K„ < 0.1, since the flow 
can still be considered as continuous. 

2.1 Boundary Conditions 

For the slip flow, the continuous models can be used if the no-slip boundary condition is modified. 
A few models have been suggested to estimate the nonzero velocity at a boundary surface^^'-*^^. 
In this study, we adopt the approach based on Taylor's expansion of the velocity around the 
wall. Thus, the first order approximation yields u \wall= Kndu/dy (positive for the inner normal 
as y = n). Consequently, the mean (basic) velocity profile is given by 

u=l-y^ + 2Kn (4) 

for —1 < y < 1. Boundary conditions for or D(f) arc not the same as previous no-slip approach, 
i.e., ^(±1) = D(p{±l) = and shall be introduced below. 
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2.2 Interface Treatment 

We consider a two-dimensional layer (slab) of uniform thickness filled with a homogeneous normal 
fluid (Newtonian viscous fluid; its dynamics is described by Navier-Stokes equations). The upper 
and lower boundaries of the layer are superfluids which are rather flexible, on which are imposed 
travelling sinusoidal waves of small amplitude a (due to SAW or peristaltic waves) . The vertical 
displacements of the upper and lower interfaces {y = h and —h) are thus presumed to be rj and 
—77, respectively, where rj = a cos ^{x — ct), X is the wave length, and c the wave speed, x and y 
are Cartesian coordinates, with x measured in the direction of wave propagation and y measured 
in the direction normal to the mean position of the interfaces. It would be expedient to simplify 
these equations by introducing dimensionless variables. We have a characteristic velocity c and 
three characteristic lengths a, A, and h. The following variables based on c and h could thus be 
introduced : 

, x , y f u , V f rj f t/j ct , p 

a? = T, 2/ = T, u = -, V = -, ri = -, = —, t = —, p = 

n n c en ch h p& 

where is the dimensional stream function. The amplitude ratio e, the wave number a, and 
the Reynolds number Re^ are defined by 

a 27r/i ch 

e=-, a = — — , i?ec = — . 
h A V 

After introducing the dimensionless variables, now, for the ease and direct representation of our 
mathematical expressions in the following, we shall drop those primes (') in those dimensionless 
variables and treat them as dimensionless. We seek a solution in the form of a series in the 
parameter e : 

^ = ^0 + e^i + e ^2 + ■ ■ ■ , ^ = (^)o + + e (^)2 + • ■ ■ , 

ox ox ox Ox 

with n = dip/dy, v = —dijj/dx. The 2D (x- and y-) momentum equations and the equation of 
continuity for the normal fluid could be in terms of the stream function if the pressure (p) 
term is eliminated. The final governing equation is 

and subscripts indicate the partial differentiation. If we presume originally the fluid is quiescent; 
this corresponds to a free pumping case and finally the velocity profile of the fluid is symmetric 
with respect to the centerline of the plane channel bounded by the superfluids. Equation above, 
together with the condition of symmetry and a uniform constant pressure-gradient in the x- 
direction, (9p/(?;r)o=constant, yield : 

V^o = i^o[(l + 2Kn)y-^], Ko = ^{-^)o, (6) 

Kq is in fact a necessary pumping to sustain a plane Poiseuille flow (pressure-driven case). Vo 
corresponds to the solution of 

d 1 

— vVo + V'oj/VVox - ^oxV^i^oy = ;^vVo, (7) 
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and 

^1 = ^{</'(y)e'"^"-*^ + </.*(y)e--(--*)}, (8) 

where the asterisk denotes the complex conjugate. The normal fluid is subjected to boundary 
conditions imposed by the symmetric motion of the wavy interfaces and the slip condition at 
interfaces. The basic slip flow now has this form (cf [16]), as u = dil^o/dy, 

« = 1 - y2 + 2Kn, (9) 
where c is the phase speed of the SAW, Kn=mfp//i. Boundary conditions become 

(f)y{±l) ± (f)yy{±l)Kn = 2Ko{l ± Kn), </.(±l) = ±1. (10) 

2.3 Numerical Approach 

The eigenvalue problem raised above could be solved by using the verified code^^ , which used the 
spectral method^'^ based on the Chebyshev-polynomial-expansion approach, once the equation 
and boundary conditions are discretized. For instance, we have, from equation (3), as a finite- 
sum approximation (reduction from oo to N), 

N 
n=0 

where Tn{z) is the Chebyshev polynomial^^'^° of degree n with z = cos{0). T„(z) are known to 
satisfy the recurrence relations 

zTn{z) = ^[Tn+liz) +Tn-liz)]. 

After substituting (p into (3) and with tremendous manipulations, we obtain the algebraic equa- 
tion 

E \P"(P" - 4)' - 3nV + 3nV - pn\n' - 4f]a,- 



p=n+4 
psn(mod2) 

N ^ ^ 

E { [2a^ + ^mi?e(4 Mof-iC-MoCn- MoCn-iMp^ - n^) - -iaReMoCnp\p^ -{n + 2f]- 

p=n+2 
p=n(mod 2) 

^iaReModn-2P\p^ - (n - 2f]}ap + iaReMon{n - l)a„ + {a'^ + iaRe[{Mof - C)a^ - 2Mo]}cnan- 

^za^i?eMo[cn_2an-2 + CniCn + Cn-l)a„ + Cnan+2] = (11) 

for n > 0, / = 1 + 2Kn, where c„ = if n > 0, and c?„ = if n < 0, (i„ = 1 if n> 0. Here, 
Mo = 1, C is the complex eigenvalue. The boundary conditions become 

N N 2/ 2 _ 1 \ 

E «n = l' E [nVKn "^"^ ' W = 2K,. (12) 

n=l n=l 

nsl(mod2) n=l(mod2) 
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The matrices thus formed are of poor condition because they are not diagonal, symmetric . 
Thus, before we perform floating-point computations to get the complex eigenvalues, we pre- 
condition these complex matrices to get less errors. Here we adapt Osborne's algorithm to 
precondition these complex matrices via rescaling, i.e., by certain diagonal similarity transfor- 
mations of the matrix (errors are in terms of the Euclidean norm of the matrix) designed to 
reduce its norm. The details of this algorithm could be traced in [18,21-22]. The form of the 
reduced matrix is upper Hessenberg. Finally we perform the stabilized LR transformations 
for these matrices to get the complex eigenvalues (please see also [18,22] for the details). The 
preliminary verified results of this numerical code had been done for the cases of Kn=0 (no-slip 
boundary conditions) in comparison with the bench-mark results of Orszag's^^ . For example, 
for Re = 10000.0, a = 1.0 of the test case : plane Poiseuille flow, we obtained the same spectra 
as 0.23752648-F i 0.00373967 for Cr+ i Cj which Orszag obtained in 1971. 

3 Results and Discussion 

After careful verification, we proceed to obtain (through tremendous searching using double- 
precision machine accuracy) the detailed spectra for the illustration of the stability of the slip flow 
in normal fluids confined between wavy interfaces. To demonstrate some of the calculated spectra 
(C,., Ci) near the regime of {Recr, a), we plot Fig. 1 by selecting 2 pairs of {Re, a)=(1562, 1.156) 
and (2982.3, 1.0783) with the corresponding Kq= 1, and 0.5 for the same Knudsen number 
(Kn=0.01). Once Ci > 0, the instability occurs! The onset of instability is easy to occur once 
the Reynolds number, wave number, or Knudsen number perturb a little again near this regime 
{Ci becomes zero and then positive). 

We then plot the neutral stability boundary curves for different cases in Fig. 2. It is clear 
that each curve is composed of two branches (one is upper and the other is lower, and they 
coalesce into a critical point {Rcc- and acr) as the Reynolds number decreases). We tune 
the Kq parameter to be 1 and 0.5, with the corresponding Knudsen number being 0.001 and 
0.01, respectively. Once the Knudsen number is set to be zero and there is no SAW effect, we 
recover the curve obtained by Orszag^^ {Recr ~ 5772). Otherwise, the resulting critical Reynolds 
numbers {Recr) are 1441,1562, 2664, 2982.3, respectively. It seems the effect of SAW propagating 
along the interface is the dominated one and will degrade the flow stability significantly. The 
slip velocity effect is minor and adverse (delay the transition). 

To understand the stability behavior related to the decay or amplification of the perturbed 
disturbance waves in the finite time for certain mode, we also illustrate their time evolution 
patterns by selecting the least unstable mode. As illustrated in Figs. 3 and 4 for Re=1441, 
a = 1.175 {Kq = 1) and Re=2664, a = 1.105 {Kq = 0.5), we can observe the oscillating or 
amplifying pattern just after a finite time (time is dimensionless and the Knudsen number is 
the same, Kn=0.001). The original disturbance (wave) will not decay for these unstable modes 
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{Cr,Ci)r^ (0.382,0.000002), and ~ (0.324,0.00000046), respectively. 

We can finally conclude that various kinds of interface noises (as illustraed here, slip velocities 
and the propagating surface acoustic wavs) will premature any instability mechanism considering 
the temporal growth of the disturbances. We have obtained more clues about the slip flow 
(which is in a non-equilibrium state) instability of quantum fluids (above their critical transition 
temperature) by considering more realistic interface conditions. Once we know the viscosities 
and/or densities of these quantum fluids, based on the obtained critical Reynolds number, we 
can then determine the critical velocity^'^ for each case. Meanwhile, these results will help 
researchers to understand the formation or generation of vorticity waves and then the route to 
low-temperature turbulence in quantum fluids. It seems the range of wave numbers relevant 
to the SAW propagating along the flexible interfaces, the Knudsen numbers and the Reynolds 
numbers of basic slip flows of normal fluids must be carefully selected for the optimal flow control 
usage in SAW applications^'^ to the investigation of ^He. Our further study will be relevant to 
those more complicated issues 2-3,22-25 
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Fig. 2 Illustration of the temporal spectra (Cr, Ci) for disturbance 

waves due to interface (i^o) and slip velocity (Kn=0.01) effects. 

Re = 1562,2982.3 for corresponding Kq = 1,0.5 and a = 1.156, 1.0783, respectively. 
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Fig. 3 Comparison of wavy interface (Kq) and slip velocity (Kn) effects 
on the neutral stability boundary of the basic flow. Kn= mfp/ft,. mfp is 
the mean free path of the quantum fluid. Rccr ~ 1441, 1562, 2664, 2982.3 
for Ko = l : Kn=0.001, 0.01, and Kq = 0.5 : Kn=0.001, 0.01. 
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Fig. 3 Illustration of the temporal evolution for disturbance 
waves due to interface (Kq) and slip velocity (Kn) effects. 

Re = 1441 for corresponding Kq=1 and q;=1.175, Kn=0.001. Time is dimensionless. 
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Fig. 4 Illustration of the temporal evolution for disturbance 
waves due to interface (Kq) and slip velocity (Kn) effects. 

Re = 2664 for corresponding Kq=0.5 and a=1.105, Kn=0.001. Time is dimensionless. 
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